Interplay between Nitrogen Dopants and Native Point Defects in Graphene 
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To understand the interaction between nitrogen dopants and native point defects in graphene, 
we have studied the energetic stability of N-doped graphene with vacancies and Stone- Wales (SW) 
defect by performing the density functional theory calculations. Our results show that N substitu- 
tion energetically prefers to occur at the carbon atoms near the defects, especially for those sites 
with larger bond shortening, indicating that the defect-induced strain plays an important role in 
the stability of N dopants in defective graphene. In the presence of monovacancy, the most stable 
position for N dopant is the pyridinelike configuration, while for other point defects studied (SW 
defect and divacancies) N prefers a site in the pentagonal ring. The effect of native point defects 
on N dopants is quite strong: While the N doping is endothermic in defect-free graphene, it be- 
comes exothermic for defective graphene. Our results imply that the native point defect and N 
dopant attract each other, i.e., cooperative effect, which means that substitutional N dopants would 
increase the probability of point defect generation and vice versa. Our findings are supported by 
recent experimental studies on the N doping of graphene. Furthermore we point out possibilities of 
aggregation of multiple N dopants near native point defects. Finally we make brief comments on 
the effect of Fe adsorption on the stability of N dopant aggregation. 

PACS numbers: 61.72.J- , 31. 15. A- , 82.45.Jn, 61.72.uf 
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I. INTRODUCTION 



Graphene is a carbon allotrope with a two dimensional 
(2D) honeycemb lattice. Since its first successful isola- 
tion in 2004jd graphene has attracted immense attention 
because of its 2D crystal lattke,with atomic thickness and 
unique electronic structures.Ercl It has opened up exciting 



opportunities for developing nanoelectronic devicesxfln 
low dimensional systems, the chemical and physical prop- 
erties of materials can be heavily affected, by the lattice 
imperfection due to the structural defects.Q The point de- 
fects in graphene can be introduced via the formation of 
vacancies or the atomic rearrangement (e.g., Stone- Wa 
defect consisting of a pentagon- heptagon (5-7) pair).E 
Vacancies in graphene may be formed in low concentra- 
tions during the growth process. Alternatively, they may 
be created intentionally by irradiating materials with 
electrons or ions or by chemical treatments. Recent ad- 
vances in microscope technologies enable us to observe 
the structural defects in graphene at an atomic resoku 
tion by the transmission electron microscope (T|EM^£rE£l 
and the scanning tunneling microscope (STM) .tlrt 3 ] On 
the other hand, recent density functional theory (DFT) 
calculations show that some of structural defects can in- 
duce localized levels-close to the Fermi level (Ep), lead=. 
ing to local chargingE-3 and/or local magnetic moments.E3 
Therefore, defects are expected to play key roles in the 
chemical functionality and electronic transport proper- 
ties of graphene-based materials. 



One simple approach to further tailor the electronic 
properties of graphene is the incorporation of het- 
eroatoms. For instance, substitution of carbon with 
nitrogen or boron atoms injects electron or hole carri- 
ers, respectively. Such carbon-based materials contain- 
ing some different elements, called carbon alloy slij, have 
been exploited recently to be a high-promising candi- 
date to replace a Pt-based-pa±alyst in the polymer elec- 
trolyte fuel cell (PEFC^E-Xe 2 ] N doping into graphene 
can be performed either directly during synthesis or by 
post-synthetic treatment .&l£| Even the controllable N 
doping has_been realized by NH3 annealing after ion 
irradiation! 3 !] or by NH3 plasma exposure.!!!! The pres- 
ence of large amount of defects in irradiated or plasma- 
treated graphene has been revealed by a pronounced D 
band in the Raman spectraElJi 3 -!! After N doping, the D 
band becomes more prominent, indicating that N dop- 
ing is likely to indu^e-^aetwieghgible amount of defects 
and bond disorders. ESEH'Ej'ES These suggest that some 
mutual effects of structural defects and N doping exist 
during the incorporation of N into graphene. Previous 
DFT calculations show that the localized edge states play 
an important role in the stability, of N substitution in 
graphene nanoribbons (GNRs)Ell'Ea and clusters.! 3 !] Thus, 
it raises a question whether the defect-induced localized 
states also have significant impacts on the N doping of 
graphene or not. Indeed, some specific configurations of 
substitutional N next to vacancies in GNRsJjiave already 
been studied theoretically in literature. L§c2l In particu- 
lar, recent workcHtEJ deals with N doping at monova- 
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cancy (MV) and divacancy (DV). In the present work, 
we study more extensively and systematically the in- 
terplay between structural defects and N doping using 
the DFT calculations with a significantly larger super- 
cell as described in detail below. Furthermore, we dis- 
cuss briefly the effect of Fe adsorption on the stability of 
N dopant arrangement, considering that the FeN^ com- 
plex in graphene may play important roles in the oxy- 
gen r^^tktLL-reaction in the carbon based fuel cell cat- 
alyst .E3tirE3o We focus on the structural and energetics 
aspects in this paper. The detailed discussion on the 
electronic structures will be given in a separate paper. 

The remainder of this paper is organized as follows. 
In Sec. [n], we introduce the computational methods for 
the calculations of N-doped graphene (N-graphene). The 
computed interaction energies between N dopants and 
structural defects are presented in Sec. III. Finally, we 
draw conclusions in Sec. |Tv| . 



II. METHOD AND COMPUTATIONAL 
DETAILS 

DFT calculations are perforated with the PWSCF code of 
the Quantum ESERESSO suiteEa in a plane- wave ultrasoft- 
pseudopotentialo approach. The exchange-correlation 
functional is treated by the generalized gradient approx=. 
imation (GGA) after Perdew, Burke, and Ernzerhof.L^I 
Spin polarization is taken into account if it exists. The 
kinetic energy cutoffs for wave function and charge are 
set to 35 and 350 Ry, respectively. A supercell con- 
structed by the 9x9 extension of the hexagonal unit cell 
of graphene with the calculated lattice constant of 2.463 
A for perfect graphene is employed to study the struc- 
tural defects and the substitution of carbon by nitrogen. 
To avoid the spurious interaction between graphene lay- 
ers, the vacuum thickness in the supercell is set to 12.0-4. 
A 3 x 3 x 1 Appoint grid in the Monkhorst-Pack schemea is 
employed to sample the Brillouin zone (BZ) of the above 
supercell. During geometry optimization, all atoms are 
fully relaxed until residual forces on constituent atoms 
are smaller than 0.01 eV/A. 

In this work, the most common types of native point 
defects (NPDs) are studied as representatives of struc- 
tural defects in graphene, namely, MV, DV, and Stone- 
Wales (SW) defect. Additionally, the hydrogenated 
monovacancy (H-MV) and the reconstructed DVs are 
also studied. The optimized atomic structures of defec- 
tive graphene are shown in Fig. [l| For DVs, besides 
the normal 5-8-5 configuration [5-8-5 DV, Fig. |l](c)], we 
consider two variants. The 555-777 configuration [555- 
777 DV, Fig. ||(d)] can be reconstructed from the 5-8-5 
DV by in-plane 90° rotation of either bond indicated by 
arrows shown in Fig. |l|(c). Similarly, the 5555-6-7777 
configuration [5555-6-7777 DV, Fig. |l](e)] can be recon- 
structed from the 555-777 DV by rotating one of the three 
bonds indicated by arrows in Fig. |l|(d). For N doping, 
we consider all the single N substitution of inequivalent 



carbons Ci in the defect region as marked in Fig. [j] in 
addition to carbon CO far from the defect for each defec- 
tive graphene. Here, such a doced N is denoted as Nci 
according to Kroger and Vinko 



III. RESULTS AND DISCUSSION 

A. Interaction between two substitutional N atoms 
in perfect graphene 

As a first step to examine the aggregation of substitu- 
tional nitrogens, the interaction energy Ai?N,N between 
two substitutional N atoms is calculated according to 



Ai?N.N = E- 



2N 



Eq — 2E 



(1) 



where Eo, E^, and E?2N are the total energies of the su- 
percells of graphene containing zero, one, and two sub- 
stitutional N atoms, respectively. Figure ||(a) shows our 
computed Ai?N,N for two doped nitrogens in the defect- 
free graphene as a function of their distance. The inter- 
action energy is generally positive (i.e., repulsive) and is 
found to be much higher for near N-N dimers than for 
distant ones. For example, as shown in Fig. ||(a), the 
interaction energy of the second nearest Nci-Nc3 pair 
is about 0.36 eV higher than that of the well-separated 
Nci-Ncio pair. Thus, such a configuration of two doped 
N atoms at the second nearest neighbors is unlikely due to 
its low stability, contrary to thepEecent proposal based on 
the STM images of N-graphene£j It is interesting to note 
that two doped N atoms at the third or seventh nearest 
neighbors (i.e., Nci-Nc4 or Nci-Ncs pairs, respectively! 
have very small interaction energy. In a recent studyCj 
the stabilities of Nci-Nc4 and Nci-Nqs pairs are ascribed 
to the low Coulomb repulsion due to the anisotropic elec- 
tron charge density redistribution induced by the N sub- 
stitution. Here we give a different interpretation to the 
trend seen in Fig. |2|(a). 

The quantity in the inset of Fig. ||(a) is obtained by 
integrating the local density of states from Ep to E-p + 
0.25 eV in the presence of a single N dopant at the CI 
site in defect-free graphene. As the impurity resonance 
state appears just above E-p as shown in the density of 
states in Fig. ^|(b), the above quantity corresponds to 
the weight distribution of the tail of the impurity state. 
Therefore, if the second N is doped at the site where 
the tail has a large weight, the two N atoms can form 
a stronger bond. Therefore we expect that C2, C4, and 
C8 sites are favorable for the second N dopant. However 
the Nci-Nc2 pair has too strong of a core-core repulsive 
interaction due to the shortest N-N distance. As the 
fourth, fifth, and sixth nearest neighbors (i.e., C5, C6, 
and C7 atoms) of Nci have less weight of impurity state 
tail, the Nci-Nc5, Nci-Nc6, and Nci-Nc7 pairs have 
lower stability than the Nci-Nc4 pair. More detailed 
discussion on the electronic structures of N-N pairs will 
be given in a separate paper. Finally we point out that 
two H adatoms are stably adsorbed on graphene in a 
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FIG. 1: (Color online) The optimized atomic structures for defective graphene with (a) MV, (b) H-MV, (c) 5-8-5 configuration 
of DV (denoted as 5-8-5 DV), (d) 555-777 configuration of DV (denoted as 555-777 DV), (e) 5555-6-7777 configuration of DV 
(denoted as 5555-6-7777 DV), and (f) SW defect. C atoms in defect region are represented by black balls and other C atoms 
are represented by gray ones. The inequivalent sites for C atom in the defect region are labeled by nonzero numbers, while C 
atom labeled with stands for the site far from the defect region. The arrows in (c) and (d) indicate possible rotated bonds 
for the reconstruction. 



wayil similar to a N-N pair treated here. The stability 
of two H adatoms was explained in terms of mesomcric 
effect, which is just another way to take into account the 
structure of the wave function at the Fermi level. e3 



B. Formation energy of native point defect in 
undoped graphene 

Before discussing the N doping of graphene with NPDs, 
we assess the stability of NPDs in graphene by calculating 
their formation energy AEq (d) as 



AE (d) = E d - E + n c /ic - "H/iH, 



(2) 



where and Eq are the total energies of a supercell of 
undoped graphene with and without a NPD. /iq and 



are the respective chemical potentials of C and H, which 
are taken as the total energy per C atom of graphene for 
C and half of the total energy of an isolated H2 molecule 
for H. nc (^h) is the number of C (H) atoms removed 
(attached) when NPD is formed. Our computed forma- 
tion energies are summarized along with those found in 
literature in Table |. For MV, H termination is found to 
reduce the formation energy by 1.948 eV/per H, indicat- 
ing that H-MV is energetically much more favorable. On 
the other hand for JDVs, we find, in agreement with pre- 
vious studies, qQOl3 that the most stable configuration 
is the reconstructed one with three pentagons and three 
heptagons (denoted as 555-777), rather than the original 
5-8-5 DV. It is also seen that the formation energy of 
SW defect is a few eV lower than that of MV and DVs, 
which is consistent with the results reported in Ref. [l4|. 
All of these validate the computational setup chosen in 
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FIG. 2: (Color online) (a) The interaction energy between two doped N atoms in defect-free graphene. The inset shows local 
density of states (0.003 e/A 3 ) integrated from Fermi level (-Ef) to _Ep + 0.25 eV for single N dopant in defect-free graphene 
and the sites considered for the N substitution. Here Ey is the Fermi level of defect-free graphene with single N dopant, (b) 
The projected density of states (PDOS) for the p orbitals of C and N atoms in defect-free graphene with single N dopant. 



the present study. 

C. Interaction between a single N dopant and a 
native point defect 

The stable positions of N substitution in defective 
graphene are explored by considering all inequivalent 
sites around the defect as shown in Fig. [j] for each type 
of NPDs. The interaction energy Ai?N,d between a N 
dopant and a NPD is defined as follows: 

Ai?N,d = (-E/N+d + Eq) — (Etsi + Ed), (3) 

where i?N+d; Eg, En, and are the total energies of 
the supercells for the N-graphene with NPDs, the perfect 
graphene, the N-graphene, and the defective graphene, 
respectively. The negative value of AE^a indicates the 
attractive interaction between N dopant and NPD. Equa- 
tion (||) can also be expressed in the following way. The 
formation energy of an N dopant in a defective graphene 
A£'d(N) and the one in a perfect graphene AE'o(N) are 
given by 

A£ d (N) = ^N+d-^d + ^c-MN, (4) 
A£ (N) = B N -Bo+Atc-MN- (5) 

Then we obtain 

AE NA = A£ d (N)-A£ (N). (6) 

The right-hand side of this equation is the difference 
in the N dopant formation energy between a defective 
graphene and a perfect one. 

The results of the calculation for A2?N,d and AEa (N) 
are shown in Fig. || for a single N dopant in various defec- 
tive graphenes as a function of the average bond length 



(dc*— c) around the substitution site (C*). A2?o(N) is es- 
timated as 0.785 eV. For each type of NPDs, the CO site 
(Fig. is separated from the NPD by several C-C bonds 
and can be an approximated site in a perfect graphene. 
Therefore, Ai?N,d for the CO site is nearly zero (at most 
about -0.3 eV for the MV case) in Fig. [|. We note the 
following facts in Fig. |J 

(1) For all the stable dopant configurations, Ai?N.d is 
negatively large implying that a N dopant and a NPD 
attract each other and that the N dopant formation en- 
ergy is significantly reduced to being even exothermic by 
the presence of NPD. 

(2) For the most stable site of N dopant in each of the 
defective graphene, the bond (C*-C) (before N doping) 
associated with the substitution site (denoted as C*) is 
the shortest and much shorter than the corresponding 
C-C bond (1.422 A) in the perfect graphene, suggesting 
that the N dopant tends to be located at the site with 
larger compressed strain in defective graphene. 

(3) In the presence of MV, the most stable site for 
the N dopant is the CI site, whose configur.ation may be 
called pyridinelike with a dangling a state Jl3 and more- 
over, evep if the dangling a state is passivated by H as in 
H-MV,li3 the N dopant still energetically prefers the CI 
site. 

(4) For DVs and SW defect, the N dopant energeti- 
cally prefers the vertex sites of the fivc-mcmbcred rings, 
the trend beingj-abserved in the N-doped carbon conju- 
gated materials, E3 which is ascribed to higher aromaticity 
according to Hiickel's rule. 

Most importantly our results suggest that the creation 
of defects in graphene before introducing N will enhance 
the incorporation of N into graphene. This would support 
the recent experimental studies, in which the authors re- 
ported that the NH3 annealing of graphene after N + -ion 
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TABLE I: Formation energy [AE (d), in eV] of MV, H-MV, DVs, and SW defect in undoped graphene. 



Defect type 


MV 


H-MV 




DV 




SW 


Configuration 


5-9 


5-9 


5-8-5 


555-777 


5555-6-7777 


55-77 


A.Eo(d) (This work) 


7.536 


5.588 


7.442 


6.616 


6.963 


4.875 


A.Eo(d) (Ref. 3) 


7.3-7.5 




7.2-7.9 


6.4-7.5 


7.0 


4.5-5.3 



0.5 
0.0 
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> -1.0 
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FIG. 3: The interaction energy (Ai?N,d, left y axis) between single N dopant (m = 1) and NPD and the formation energy 
of N substitution [AJS d (N), right y axis] in defective graphene with (a) MV, (b) H-MV, (c) 5-8-5 DV, (d) 555-777 DV, (e) 
5555-6-7777 DV, and (f) SW defect as function of average C*-C bond length (dc*-c) around the substitution site (C*) before N 
doping. The open hexagonal symbols stands for C* at a six-membered ring, the solid pentagonal ones for C* at a five-membered 
ring, and the up-triangle symbols for CO* far away from the defect region. The positions of substitution sites Ci* (here i is the 
number used to mark the inequivalent sites around the defect region) can be referred to Fig. [I]. 
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irradiatior£3 or the NH3 plasma exposures!! of graphene 
can be used to realize the N doping of graphene in a 
controllable manner. 

So far, Ai?N,d has been discussed from the view point 
of the effect of NPD on the N dopant formation energy. 
Similarly, Ai?N,d can be regarded as the difference in the 
formation energy of defects after and before N doping. 
More concretely, the NPD formation energy in the pres- 
ence of N dopant AEn (d) is given by 

A£ N (d) = E N+d - .En + ncvc - nnfJ-n- (7) 

The one in a perfect graphene AE (d) is given by Eq. (||). 
Then we obtain 

A£ Njd = AE N (d)-AEo(d). (8) 

The results depicted in Fig. || can now be regarded as 
the energy gain in the NPD formation energy caused by 
the presence of N dopant. Therefore, they suggest that N 
doping would increase the probability of point defect gen- 
eration in graphene. In experimental studies, the Raman 
spectra of N-graphene have an intense D band, which is 
generally .asmhecLto the defects in the sp 2 C network of 
graphene»yy 

D. Multiple N dopants near native point defects 

For defect-free graphene, as the interaction between 
two N dopants is repulsive [Fig. 0(a)], aggregation of N 
dopants in a small region is unlikely. However, we showed 
that a N dopant and a NPD attract each other quite 
strongly. Therefore, we expect that NPD may induce 
N dopant aggregation. Figure ^ presents the most stable 
configurations of multiple N substitutions (m = 2,3, and 
4) in defective graphene. First we point out some char- 
acteristic features seen in the stable configurations. In 
the presence of MV, one of two N dopants occupies the 
site next to the vacancy to form a pyridinclikc N and the 
other one substitutes for a C atom at a five-membered 
ring [see Fig. ^(b)]. This configuration is stable by 0.257 
eV compared to that of two pyridinelike N atoms at a 
MV. For three N dopants at MV, all of them form the 
pyridinelike configuration. Such pyridinelike N at MV 
in carbon nanotube (CNT) is proposed to be responsi- 
ble for the introduction of a large electron donor state 
in N-doped CNT.Ej The present study shows that the 
configuration of the three pyridinelike N atoms at MV 
has a very high stability against the others. For exam- 
ple, the total energy difference between the most stable 
configuration [see Fig. @(c)] and the second stable one 
[Nci + Nc7 + Nc8, where C8 is the second nearest neigh- 
bor of CI as shown in Fig.|](c).] of three N substitutions 
at MV is about 1.763 eV. For two N dopants at the 5-8- 
5 DV the most stable configuration corresponds to two 
pyridinelike N atoms [see Fig. ||(f)]. This configuration 
is more stable than the one with two N dopants at C4 
sites by 0.223 eV. These two pyridinelike nitrogens tend 



to be closer together, contrary to the corresponding ni- 
trogens at the armchair edge of GNR.E3 Up to three N 
dopants at the 555-777 DV or up to two N dopants at the 
SW defect, the most stable configuration is the combined 
single N dopants, where substitutional N atoms occupy 
the vertex sites in different five-membered rings. 

Now we discuss the possibility of N dopant aggrega- 
tion near a NPD in a quantitative way. We generalize 
Eq. (§) to the case of multiple N dopants where m N 
atoms substitute for m C atoms near a NPD to obtain 

A£ l d (mN) = E mN+d — Ed + mp,c — to/zn- (9) 

Figure ||(a) presents AEd(mN) for the configurations in 
Fig. H It is clear that N dopants are particularly stable 
at MV and that four N dopants at 5-8-5 DV shows also 
strong stability. The formation energy of a complex of 
mN dopants plus a NPD is given by the sum of the NPD 
formation energy AE (d) of Table | and A£ d (mN). Al- 
though the formation energy of a NPD without N doping 
is smallest (largest) for SW (MV), 3N doping at the NPD 
makes MV have the smallest formation energy. Further- 
more, the formation energy of 5-8-5 DV, which is the 
second largest in undoped graphene, sharply decreases 
with N doping and becomes the second smallest with 4N 
doping. 

In addition to these quantities, we define the formation 
energy of ith N dopant in a system containing already 
(i — 1) N dopants near a NPD by the following equation: 

Ai?d(Nj) = £iN+d - E(i-i)N+d + nc - Mn- (10) 

Clearly, AE d (Ni) is the finite difference of AE d (iN) with 
respect to i and the following relation holds. 

m 

A£ d (mN) = ^A£ d (N i ). (11) 

i=l 

Figure |(b) shows AE d (N t ) for i ranging from one to four 
for each of the different NPDs. So long as AE d (Nj) is 
smaller than 0.785 eV, which is the formation energy of 
N dopant in perfect graphene indicated with the dotted 
line in Fig. |^(b), an aggregate of i N dopants at the NPD 
is energetically stable if other NPDs are not present. 

However, the situation may be different in the pres- 
ence of other NPDs. We consider the case where two 
NPDs, designated as di and d2, exist before N doping 
and calculate the following energy 

£7(ftiN + di, fe 2 N + d 2 ) - AE dl (kiN) + AE d2 (fc 2 N), 

which is the sum of formation energies of k\N dopants at 
a di-NPD and £; 2 N dopants at a d 2 -NPD. We then find k\ 
(/c 2 ) which minimizes the above energy for a given set of 
di , d 2 and m = k\ + fc 2 . We only present some examples 
of this analysis in the following. 

First, we consider the following cases of equiva- 
lent two NPDs and present the calculated values of 
J5(jfeiN + di, fc 2 N + d 2 ) in Table Q. 
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FIG. 4: (Color online) The optimized structures for the most stable configurations of m N dopants (m = 1, 2, 3, and 4) in 
defective graphene with MV [(a)-(d)], 5-8-5 DV [(e)-(h)], 555-777 DV [(i)-(l)], and SW defect [(m)-(p)]. The blue and light 
gray balls stand for the N and C atoms, respectively. 



TABLE II: The sum of formation energies [_E(fciN + di, feN + d2), in eV] of fciN dopants at a di-NPD and feN dopants at a 
d2-NPD for the case of di = d2. 



di = d2 m 


= fci + fc 2 


£(fc 1 N + di,faN + 


d a ) 






ki = 1 


fci =2 


fci = 3 


fci = 4 


MV 


2 


-4.298 


-2.865 






MV 


3 


-5.015 


-5.015 


-4.136 




MV 


4 


-6.285 


-5.731 


-6.285 


-4.373 


5-8-5 DV 


2 


-1.012 


-1.467 






5-8-5 DV 


3 


-1.973 


-1.973 


-2.065 




5-8-5 DV 


4 


-2.572 


-2.934 


-2.572 


-3.657 


555-777 DV 


2 


-1.035 


-0.924 






555-777 DV 


3 


-1.442 


-1.442 


-0.818 




555-777 DV 


4 


-1.336 


-1.849 


-1.336 


0.008 


SW 


2 


-0.237 


-0.535 






SW 


3 


-0.654 


-0.654 


-0.101 




SW 


4 


-0.220 


-1.071 


-0.220 


0.993 
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FIG. 5: (Color online) (a) The formation energy of N dopants [AEd(mN)] in graphene with a NPD of some different kinds as 
a function of the number of substitutional N atoms, (b) The formation energy of ith N dopant in a system containing already 
(i — 1) N dopants near a NPD [A£d(Ni)] for i in the range from 1 to 4. The formation energy of single N dopant in perfect 
graphene is indicated with the dotted line in (b). 



(1) di = d 2 = SW: For m = 2, (fei = 2, k 2 = 0) or 
equivalently (fci = 0, k 2 = 2) is more stable than (ki = 1, 
k'2 = 1). For m=3, (fci = 2, k 2 = 1) or equivalently 
(k\ = 1, ki = 2) is more stable than (k\ = 3, = 0) and 
(fci = 0, k 2 = 3). Therefore aggregation of 2 N dopants 
is possible. 

(2) di = d 2 = 5-8-5 DV: For to — 2, the situation is the 
same as in SW. For to = 3, the situation is opposite to the 
one in SW. For m = 4, (ki = 4, k2 — 0) or equivalently 
(ki = 0, &2 = 4) is more stable than other partitions 
of 4 N dopants. From this analysis, we conclude that 
aggregation of 4 N dopants at one 5-8-5 DV is possible 
even if more 5-8-5 DVs may exist before N doping. 

(3) di = d 2 = MV: As the formation energy of the 
first N doping is strongly negative, for any m > 1, (k\ = 
1, &2 = m — 1) is not stable among other partitions. 
Therefore, if many MV may exist, the configuration with 
the maximum number of MVs with single N dopant will 
be most stable. 

For inequivalent two NPDs (see Table SI in the Sup- 
plemental Material at the link in |5^), if one of the two 
NPDs (di, for example) is MV, selective N doping will 
occur at MV due to the large negative formation energy 
of N dopants at MV. However, in our analysis, if the other 
NPD is 5-8-5 DV or 555-777 DV and m = 4, (fci = 3, 
&2 = 1) is slightly more stable than (fci =4, k 2 = 0). 

Finally we discuss the effect of the adsorbed Fe atom 
on the stabilities of the m pyridinelike N dopants at the 
MV [m — 1 and 3, as shown in Fig. [|(a) and (c)] and 
5-8-5 DV [m = 2 and 4, as shown in Fig. ||(f) and (h)]. 
Table HI is a summary of calculated results. Though de- 



tailed discussion on the physical and chemical properties 
of Fe adsorption to N-doped graphene will be given sep- 
arately, we here point out an important consequence of 
Fe adsorption on the stability of N dopant configuration. 
Combining AS (d) of Table §, AE^(mN) in Fig. |(a) aad 
E a d (Fe) in Table [II , we find that the formation energyLD 



of four N dopants at 5-8-5 DV becomes lower than that 
of three N dopants at the MV by 2.357 eV in the presence 
of adsorbed Fe, while the former is higher than the latter 
by 0.385 eV in the absence of adsorbed Fe. This suggests 
that four pyridinelike N dopants at the 5-8-5 DV may be 
quite stable in the presence of iron. The strong stability 
of the FeN* complex was pointed out also in the carbon 
nanotubeM3 



IV. CONCLUSION 

We have studied the effect of native point defects on 
the N doping of graphene using DFT electronic structure 
calculations. Our calculations show that substitutional 
N doping tends to occur at the carbon sites with larger 
shrinkage of the defect-induced bond and also at the ver- 
tex site of a pentagonal ring if it exists. The presence of 
native point defects can lower the formation energy of N 
dopant from being endothermic in defect-free graphene 
to being exothermic in defective graphene. This suggests 
that the intentional creation of defect before introducing 
N dopant will enhance the N doping of graphene. On the 
other hand, the formation energy of defects is reduced 
after N doping, indicating that N doping would increase 
the probability of point defect generation in graphene. 
We also analyzed possibilities of multiple N doping at 
a NPD due to the strong attractive interaction between 
a native point defect and a N dopant. In the actual N 
doping process of graphene, the partial pressure of N2 or 
NH3 gas, the high-temperature heat treatment, and the 
kinetic factors may also affect the distribution of doped 
N. Nevertheless, our analysis may give some insight into 
the stability of N dopant configuration in the defective 
graphene. In the present paper, we focused on structural 
and energetics aspects of the problem. The detailed dis- 
cussion on the underlying electronic structures will be 
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TABLE III: FeNz centers are formed by the Fe adsorbate and the pyridinelike N at MV and 5-8-5 DV in graphene: adsorption 
energy [_E a( j(Fe), in eV/atom] of Fe atom, total magnetic moment (M tot , in /iB), local magnetic moment of Fe atom (rriFe, 
in /iB/atom), bond length between Fe and its nearest neighbors [c?Fe-N or dFe-c, in A], and height of Fe atom (zf c , in A) 
with respect to the atomic plane of graphene sheet. The adsorption energy of Fe atom is calculated as follows: _E at j(Fe) = 
_Ep e+m N + d — Ed — /iFc, where i?Fe+mN+d is the total energy of FeN^ center embedded in graphene. /iFe is the chemical potential 
potential of iron and here it is taken as the total energy of a Fe atom. 



Configuration 


£ ad (Fe) 


M tot 


m-pc 




dFo-C 


ZFc 


IN at MV 


-5.425 


0.96 


0.83 


1.744 


1.763 


1.574 


3N at MV 


-4.837 


3.40 


3.11 


1.880 




1.789 


2N at 5-8-5 DV 


-7.230 


2.30 


2.70 


1.959 


1.879 


0.0 


4N at 5-8-5 DV 


-7.579 


1.98 


2.06 


1.895 




0.0 



given in a future publication. 
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